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Abstract 

For the experimental search of neutralino dark matter, it is important to know its allowed mass 
and scattering cross section with the nucleon. In order to figure out how light a neutralino dark 
matter can be predicted in low energy supersymmetry, we scan over the parameter space of the 
NMSSM (next-to-minimal supersymmetric model), assuming all the relevant soft mass parameters 
to be below TeV scale. We find that in the parameter space allowed by current experiments the 
neutralino dark matter can be as light as a few GeV and its scattering rate off the nucleon can reach 
the sensitivity of XENON100 and CoGeNT. As a result, a sizable parameter space is excluded by 
the current XENON100 and CoGeNT data (the plausible CoGeNT dark matter signal can also 
be explained). The future 6000 kg-days exposure of XENON100 will further explore (but cannot 
completely cover) the remained parameter space. Moreover, we find that in such a light dark 
matter scenario a light CP-even or CP-odd Higgs boson must be present to satisfy the measured 
dark matter relic density. Consequently, the SM-like Higgs boson hsm may decay predominantly 
into a pair of light Higgs bosons or a pair of neutralinos so that the conventional decays like 


hsm > yy is much suppressed. 


PACS numbers: 12.60.Jv,11.30.Qc,12.60.Fr,14.80.Cp 


Introduction: Experiments for the underground direct detection of cold dark matter X 
have recently made significant progress. While the null observation of x in the CDMS 
and XENON100 experiments has set rather tight upper limits on the spin-independent (SI) 
cross section of X-nucleon scattering Al d, the CoGeNT experiment [3] reported an excess 
which cannot be explained by any known background sources but seems to be consistent 
with the signal of a light Y with mass around 10 GeV and scattering rate around 1074 
cm?. Intriguingly, this range of mass and scattering rate is compatible with the dark matter 
explanation for both the DAMA/LIBRA data and the preliminary CRESST data [4]. So 
far, due to the inconsistency of the CoGeNT result with the CDMS or XENON result, it 
is premature to draw any definite conclusion about the existence or nonexistence of a light 
xy. However, considering much effort is being paid on the search of a light x in experiments, 
it is theoretically important to check the possible new physics prediction for a light x and 
examine its related phenomenology (such as the Higgs boson search) at the LHC. In this 
work we will focus on low energy supersymmetry, where the lightest neutralino X? serves as 
the dark matter candidate, and perform an intensive study of the light ¥? scenario. 

The most popular model for low energy supersymmetry is the minimal supersymmetric 
standard model (MSSM). In this model, we find from our scan that the neutralino ¥? must 
be heavier than about 28 GeV. The main reason for the absence of a lighter X? is its 
dominant annihilation channel is X9x9 — bb through s-channel exchange of the pseudoscalar 
Higgs boson (A) and the measured dark matter relic density requires m4 ~ (90-100) GeV 
and tan@ ~ 50, which is in conflict with the constraints from the LEP and B physics 
experiments ld Here we emphasize that for the above region the effects of the charged 
Higgs on B > X,y are unacceptably large, and even with a fine tuning of the contributions 
from the stop/chargino diagrams, such large effects cannot be reduced to_an acceptable 
level. Our results for the MSSM are in agreement with Id. but differ from [8] in which the 
considered constraints, such as the invisible Z-decay and the productions of neutralinos or 
Higgs bosons at LEP II, are weaker than in our study. Our conclusion also differs from 
because we used more accurate formula in calculating the process B > Xyy 10 . Since 
the neutralino dark matter in the MSSM cannot be so light as suggested by the CoGeNT 
data (albeit not corroborated by XENON100 or CDMS), here we do not present our MSSM 
results in detail. 


Another popular model for low energy supersymmetry is the next-to-minimal supersym- 


metric standard model (NMSSM) which extends the MSSM by adding one gauge singlet 
chiral superfield S. This model is well motivated because it provides a solution to the u- 
problem and the little hierarchy problem the MSSM suffers from. For this model we perform 
an intensive scan over its parameter space by assuming all the relevant soft mass parameters 
below TeV scale and considering various experimental constraints. We find that in this model 
the neutralino dark matter can be as light as a few GeV and its spin-independent scattering 
cross section with the nucleon can reach the sensitivity of CoGeNT and XENON100. 

We emphasize that our study is not restricted to explain any special experimental results 
like the CoGeNT or DAMA/LIBRA data. Instead, we aim to investigate the characteristics 
of a light ¥?, such as its lower mass bound and scattering rate off the nucleon, and also 


examine the related Higgs phenomenology at the LHC. 


The NMSSM: We start our analysis by recapitulating some basics of the NMSSM. Its 
superpotential and the associated soft-breaking terms in the Higgs sector are given by 


AA A L a 
W = AĤ, Ha + 345", (1) 
Vsoft = mig Hal? + mg, ul? + mg] SÅ 


+ (AA Hy Ha + hic.) + (74.5 + he), (2) 


where Hy, H,, and S denote scalar components of the superfield Ha, H, and $, respectively. 
After using the minimization condition of the Higgs potential, this sector is described by 
three dimensionless parameters (tan 8, Å, x) and three dimensionful parameters (u, A), 
A,,). Due to the imposed Z3 symmetry, the superpotential does not contain dimensionful 
parameters and thus all dimensionful parameters are generated by the soft-breaking masses 
which, as required by the electroweak symmetry breaking, should be naturally below TeV 
scale hl Other free parameters are the same as in the MSSM, i.e., the soft masses for 
sfermions and gauginos as well as the trilinear soft couplings. 

Due to the presence of s , the NMSSM predicts five neutralinos, three CP-even Higgs 
bosons (hj,2,3) and two CP-odd Higgs bosons (a12) [LI]. In general, the neutralino mass 
eigenstates are the mixture of the MSSM neutralino fields and the singlino field which is 
the fermion component of s ; the CP-even (odd) Higgs mass eigenstates are similarly the 
mixture of the CP-even (odd) MSSM Higgs fields and the real (imaginary) part of the scalar 
component of $. An important feature of the NMSSM is that one of the CP-even (odd) 
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Higgs bosons may be singlet-like and thus can be very light al. This feature is particularly 

useful for light Y? scenario since it opens up new important annihilation channels for X$, 

i.e., either iim in ae of hj (or a) or into a pair of fermions via s-channel exchange of 
3 


hy (or az) Hlad 


latter case, due to the very weak couplings of hı (a1) with ¥? and with the SM fermions, a 


. For the former case, X? must be heavier than hı (a1); while for the 


resonance enhancement (i.e. mn, or Ma, must be close to 2m,o) is needed to accelerate the 
annihilation. So a light ? should be necessarily accompanied by a light hj or a, to provide 
the required dark matter relic density. 

Now we discuss how to get a light hj or a, in the NMSSM. A light a, can be easily 
obtained when the theory approaches to the U(1)z or U(1)pg symmetry limit, which can be 
realized by setting the product KA, to be negatively small fol In contrast, a light hj can 
not be obtained so easily, but, as shown below, it can still be achieved by somewhat subtle 
cancellation via tuning the value of A,. We note that for any theory with multiple Higgs 
fields, the existence of a massless Higgs boson implies the vanishing of the determinant of its 
squared mass matrix and vice versa. For the NMSSM, at tree level the parameter Å, only 
enters the mass term of the singlet Higgs bosons and thus the determinant (Det M?) of the 
mass matrix for the CP-even Higgs bosons depends on A, linearly [11]. When other relevant 
parameters are fixed, one can then obtain a light hı by varying A, around the value Å, which 
is the solution to the equation Det M? = 0. In practice, one must include the important 
radiative corrections to the Higgs mass matrix, which will complicate the dependence of M? 
on A,. However, we checked that the linear dependence is approximately maintained by 
choosing the other relevant parameters at the SUSY scale, and one can solve the equation 


iteratively to get the solution Åp. 


Numerical scan and results: In order to study light X? scenario in the NMSSM, we scan 
randomly over the parameters in the neutralino and Higgs sectors by requiring 0 < A,K < 
0.7, I < tan < 60,0 < u, Ay, Ma < 1 TeV and 0 < Mı < 100 GeV. Here the ranges of À, 
k and tan @ are determined by the perturbativity of the theory , the ranges of u and A) 
are suggested by the electroweak symmetry breaking, and the narrow range of bino mass 
M; is chosen to facilitate a light ¥?. Since the gluino mass and the soft parameters in the 
squark sector affect little on the properties of %9, we set all of them to be 1 TeV. As for 


the soft slepton parameters, since they influence the muon anomalous magnetic moment a, 


which in turn can limit the important parameter tan 8, we assign them a common scale m; 
and vary it below TeV scale. 

Since a light ¥? is very likely to be accompanied by a light hı or a,, we perform two 
independent scans aiming at a light hj and a light a, respectively. For the light hj case, we 
vary A, around A, which is obtained by solving the equation Det M? = 0; while for the 
light a, case, we simply vary Å, in the range [—200, 0] GeV. 

In our scans, we require all dimensionful parameters in the Higgs potential like my, 
and my, below TeV scale and keep the parameter points which yield myo < 20 GeV. The 
constraints considered in our scan are the following mt (1) We require X2 to account for the 
dark matter relic density 0.105 < Qh? < 0.119; (2) We require the NMSSM contribution to 
explain the deviation of the muon a,, i.e., a? — a3M = (25.54 8.0) x 1071", at 20 level; (3) 
The LEP-I bound on the invisible Z-decay, r(Z — ¥?X?) < 1.76 MeV, and the LEP-II upper 
bound on o(ete~ — X°X°), which is 5 x 107? pb for i > 1, as well as the lower mass bounds 
on sparticles from direct searches at LEP and the Tevatron; (4) The constraints from the 


direct search for Higgs bosons at LEP-II, including the decay modes h > hıhı, aja > 4f, 


which limit all possible channels for the production of the Higgs bosons; (5) The constraints 


from B physics observables such as B > X,7, Bs > pty, Ba > X,utu”, BY > ttv, 
YT — ya, the a mixing and the mass difference AM, and AM,; (6) The constraints 
from the precision electroweak observables such as piept, sin? gr mw and R,; (7) The 
constraints from the decay T(1S) > yhi > yimtn, KK), Y(nS) > yhi > yutu 
(n = 1,2,3) and the Tevatron search for a light Higgs boson via 4u and 2427 signals 
The constraints (1-5) have been encoded in the package NMSSMTools mt We use this 
package in our calculation and extend it by adding the constraints (6, 7). As pointed out in 
ha the constraint (7) is important for a light Higgs boson. 

In Fig. [we display the surviving parameter samples, which are simultaneously projected 


on the alm, plane in the left frame and the olLmn, or g’! 


-Mq, plane in the right frame. 
In our calculation we use the formula in for the scattering rate and choose fr = 0.023, 
fi) = 0.034, g = 0.019, = and fi = = 0.020 as input. Note that, 
motivated by recent lattice simulation , we take a very small value of fr, and consequently 
our estimation of the rate is rather conservative. 


The left frame of Fig. [I]clearly shows that in the NMSSM the neutralino dark matter can 
be as light as several GeV, while the right frame of Fig. [I shows that such a light neutralino 
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FIG. 1: The scatter plots of the parameter samples which survive all constraints. The samples 
denoted by ‘x’ (red), which are simultaneously projected on the om, plane in the left frame and 
the o°!-mp,, plane in the right frame, are characterized by a light hj; while the samples denoted 
by ‘A’ (green), which are simultaneously projected on the alm, plane in the left frame and the 
0”-mq, plane in the right frame, are characterized by a light a,. The curves are the limits from 
CoGeNT i CDMS [1] and XENON100 J while the contour is the CoGeNT-favored region 


The future XENON100 (6000 kg-days exposure) sensitivity is also plotted [17]. 


is accompanied by either a light hj or a light a1. So the surviving samples were classified into 
two sets, characterized respectively by a light hı and a light aı. In both the light-h, case and 
the light-a; case the surviving samples give a bino-dominant ¥? and correspondingly a small 
M; (<30 GeV). The allowed regions for other parameters are listed in Table [| We see that 
the value of u is not so large (below 300 GeV in light-h; case and 350 GeV in light-a, case) 
because a low u can enhance the coupling of x1 to hj and a, to get the correct relic density. 
Also we see a moderately loose bound on Ma, (we checked that among the surviving samples 
about 2% have Mma, > 60 GeV). This is due to a possibly large tan 8, which can enhance the 
couplings of a; to the SM fermions so that ma, may deviate from 2m o significantly. In the 
light-h, case, about 75% of the surviving samples are found to satisfy mp, < myo, in which 


the annihilation mode ¥?X¥? — hıhı plays a crucial role in getting the required dark matter 


TABLE I: The allowed parameter ranges for the light-h, case (first row) and light-a; case (second 


row). The dimensionful parameters pi, Å; and mp, q, are in unit of GeV. 


relic density. In contrast, in the light-a, case most surviving samples are found to satisfy 
Ma, 2 Mo, in which the dominant annihilation channel of dark matter is VE = at — ff. 
From Fig. [I] we see that the scattering rate of the light dark matter can reach the sen- 
sitivity of XENON100 and, consequently, a sizable parameter space is excluded by the 
XENON100 (2011) data [2]. The future XENON100 experiment (6000 kg-days exposure) 
17] can further explore (but cannot completely cover) the remained parameter space. Note 
that in the light-hj case the scattering rate can be large enough to reach the sensitivity of 
CoGeNT and can cover the CoGeNT-favored region. The underlying reason is that the x- 
nucleon scattering can proceed through the t-channel exchange of the CP-even Higgs bosons, 
which can be enhanced by a factor 1/ må, for a light hi ha while a light a, can not give 
such an enhancement because the CP-odd Higgs bosons do not contribute to the scattering 
claimed that the NMSSM is unable to 


explain the CoGeNT data because they did not consider the light-h, case. 


in this way. We noticed that the studies in 


Note that the light-h; samples are separated into three regions of mn, as shown in the 


right frame of Fig. This is due to the combined constraints from the dark matter relic 


density and the processes Y(1S) —> yhı (with hj 3 ntr, KK, pty), Ba > Xsut u 
and pp > hsm > 4p(or 2427) at the Tevatron, which tightly limited the couplings and 
mass of hj. For example, since in this case hj mainly acts as the product of dark matter 
annihilation, its coupling to Y? must be moderately large to get the correct relic density, 
while its couplings to the SM fermions and hsm must be small to suppress the rates of the 
processes mentioned above [15]. The situation is different for the light-a,; case, where the 
relic density and the LEP search for Higgs boson require ma, > 8GeV. For such a ‘heavy’ ag 
the above mentioned low energy processes give no stringent constraints and thus the light-a; 
samples are not separated into different regions of m,,, as shown in the right frame of Fig. 


Implication on Higgs physics: In the NMSSM the light x? scenario may predict rather 


peculiar Higgs phenomenology due to the presence of the light particles. Among the pre- 
dicted Higgs bosons, the SM-like Higgs boson Asm, defined as the CP-even Higgs boson with 
largest couplings to Z° pair, will be the most important one to be searched at the LHC 
since it mainly responsible for electroweak symmetry breaking. So we focus on hsm in our 


following discussion. 
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FIG. 2: Same as FigIT] but showing the decay branching ratios of the SM-like Higgs boson hgm. 
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Here Br(hsm > XlX?) denotes the total rates for all possible hsm > X¥°X2 decays. 
iXj iAj 


In the light ¥? scenario, hsm may decay exotically into L, hıhı or aja, and conse- 


quently the conventional decays are reduced. This feature is illustrated in Fig. P| which 
shows that the sum of the exotic decay branching ratios may exceed 50% and the the tradi- 
tional decays hsm > bb, rT, WW*, yy can be severely suppressed. Numerically, we find that 
the branching ratio of hsm > bb is suppressed to be below 30% for all the surviving samples 
in the light-h; case and for about 96% of the surviving samples in the light-a; case. For 
the remaining 4% samples in the light-a,; case, due to the kinematical forbiddance of the 
decay hsm > 141, the ratio of hsm — bb usually exceed 30% and may even approach its 
SM value (~ 70%). The samples with the ratio exceeding 65% are found to be characterized 
by ma, > 58 GeV, ma, > 350 GeV and tan 8 > 37. 

Another interesting feature shown in Fig. Blis that, due to the open-up of the exotic 
decays, hsm may be significantly lighter than the LEP bound. This situation is favored by 


the fit of the precision electroweak data and is of great theoretical interest [21]. 
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FIG. 3: Same as Fig. [1] but showing the diphoton production rate of the SM-like Higgs boson at 
the LHC. 


Since the conventional decay modes of hsm may be greatly suppressed, especially in the 
light-h; case which can give a rather large y — N scattering rate, the LHC search for hsm 
via the traditional channels may become difficult. Noting hsm is bounded from above by 


about 130 GeV and hence its most important discovering channel at the LHC is the di- 


photon signal, we show the di-photon rate at the LHC with vs = 14 TeV in Fig. In 
calculating this rate, we used the narrow width approximation and only considered the 
leading contributions to pp > hsm from top quark, bottom quark and the squark loops. 

Fig. B] indicates that, compared with the SM prediction, the NMSSM rate in the light 
X2 scenario is suppressed to be less than 20 fb for the light-h; case, and for the light-a; 
case most samples (about 96%) give the same conclusion. Since in the light-h, case the 
x — N scattering rate can reach the CoGeNT sensitivity, this means that in the framework 
of NMSSM the CoGeNT search for the light dark matter will be correlated with the LHC 
search for the Higgs boson via the di-photon channel. 

Note that, as shown in the left frame of Fig] a few samples can give a di-photon 
rate which is comparable with or even exceeds its SM prediction. We checked that these 
samples predict approximately same decay branching ratios of hsm as the SM Higgs boson, 
and the excess is mainly due to the slight suppression of the width of hsm — bb so that 
Br(hsm — yy) is enhanced. We also checked that the null result of the future XENON 
6000 kg-days exposure will imply a Higgs di-photon signal below 20 fb at the LHC with 99% 
probability. 

Finally, we point out that our light-h, scenario is different from the scenario considered 
in Ri The basic ideas of are: (1) Consider a special part in the parameter space, 
which is characterized by À, « > 0 and A, œ utan B so that Det M? ~ 0 to get a light 
hı; (2) Consider a nearly decoupled $ so that the singlino serves as the dark matter with 
its annihilation and its scattering with the nucleon proceeded mainly by exchanging a light 
singlet a; and hı respectively. Such a treatment obviously has the unnaturalness problem 
in electroweak symmetry breaking since the condition A, ~ utan 9 usually pushes the soft 
mass my, in Eq. up to several TeV. We note the results of also suffer from this 
problem. In our scenario, however, we keep the naturalness by requiring all soft masses to 
be below TeV scale. 

Conclusion: We scrutinized the light neutralino dark matter scenario in the NMSSM 
by scanning over the parameter space with all the relevant soft masses below TeV scale. 
We found that in the parameter space allowed by current experiments the neutralino dark 
matter can be as light as a few GeV and its scattering rate with the nucleon can reach 
the sensitivity of XENON100 and CoGeNT (the CoGeNT signal can be explained). The 
present XENON100 and CoGeNT data can exclude a large parameter space, and the future 
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6000 kg-days exposure of XENON100 can further explore (but cannot completely cover) the 
remained parameter space. In such a light dark matter scenario, a light CP-even or CP-odd 
Higgs boson must be present to satisfy the measured dark matter relic density. As a result, 
the SM-like Higgs boson Asm may dominantly decay into a pair of light Higgs bosons or a 
pair of neutralinos, and consequently the conventional decays like the di-photon signal at 


the LHC will be much suppressed. 
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